Quantum chaos and effective thermalization 
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We demonstrate eftective equilibration for unitary quantum dynamics under conditions of classical 
chaos. Focusing on the paradigmatic example of the Dicke model, we show how a constructive 
description of the thermalization process is facilitated by the Glauber Q or Husimi function, for 
which the evolution equation turns out to be of Fokker-Planck type. The equation describes a 
competition of classical drift and quantum diffusion in contractive and expansive directions. By this 
mechanism the system follows a 'quantum smoothened' approach to equilibrium, which avoids the 
notorious singularities inherent to classical chaotic flows. 



How do thermally isolated many particle quantum sys- 
tems relax towards stationary states? Along with the 
recent advances in quantum optics and cold atom ex- 
perimentation, this fundamental question of theoretical 
physics presents itself under new perspectives. The inter- 
acting, and in many instances thermally decoupled many 
particle states realizable in cold atom devices show novel 
forms of relaxation which can be monitored at hitherto 
unimaginable degrees of resolution. Examples of uncon- 
ventional relaxation phenomena include the formation of 
negative temperature states [I], the buildup of textures in 
Brillouin zones by an interplay of interaction and single 
particle kinematics f5^, or the relaxation into ergodic and 
other types of 'non-thermal' distributions [3HS] . 

There appears to be a concensus that for fully devel- 
oped chaos, and in the absence of conservation laws, 're- 
laxation' towards equipartition on a shell of conserved 
energy ('microcanonical distribution') represents a uni- 
versal paradigm. Recent work describes the ensuing equi- 
librated state in terms of energy eigenfunctions whose 
phase space profiles are presumed ergodic [SHS]. While 
computationally efficient, this scheme does not reveal the 
physical processes leading to thermalization. 

In this letter we approach the phenomenon from a 
more dynamical perspective. We present our case for 
the paradigmatic example of the Dicke model [6 which 
describes the coupling of a (large) spin and an oscilla- 
tor. At a critical coupling strength, the model undergoes 
a quantum phase transition into a 'superradiant' phase, 
characterized by a nonvanishing mean oscillator ampli- 
tude and chaotic dynamics [7 . In view of the recent ex- 
perimental observation of that phase transition j^, the 
issue of thermalization has become concrete. 

Our approach hinges on a manageable description of 
the impact of quantum fluctuations on the semiclassical 
dynamics. A representation of the quantum density op- 
erator, by Glauber's Q-function (aka Husimi function), 
turns out to be the key. Von Neumann's equation for the 
density operator assumes the form of a Fokker Planck 
equation for Q. The equilibration processes can be then 
understood by appreciating the competition of classical 
drift and quantum diffusion in expansive and contractive 



directions of the 'flow' supporting Q. We augment the 
general discussion by an exact solution of a toy model 
which arguably reflects the essence of the full problem. 

In the end we will indicate that our methodology is 
applicable to a whole class of chaotic dynamics. Even 
certain periodically kicked systems qualify. 

Dicke Model — . We write the Dicke Hamiltonian as 

H = h^LOoJ;, +uja''a + g.^{a + a^)X^ . (1) 

Here, Ja,a = x,y,z are spin operators acting in a 
spin-j representation, and a/a^ are photon annihila- 
tion/creation operators. The first two terms in ([T]) re- 
spectively describe spin precession about the Jz-axis with 
frequency wq and harmonic oscillation with frequency w. 
The last term describes the coupling of spin and oscillator 
with coupling constant g. Crucially, the interaction con- 
tains the so called antiresonant terms J+a^ + J-a where 
Jj. — Jx±iJy are the familiar raising and lowering oper- 
ators. This fact makes the model distinct from an inte- 
grable variant where these terms are neglected (rotating 
wave approximation). At a critical value of the coupling 
strength, gc = y/uJOJo/2, the non-integrable model ([T]) 
undergoes a quantum phase transition into a 'superradi- 
ant phase'. For values g > g^ the photon 'amplitude' a 
builds up a nonvanishing excitation value, and the dy- 
namics becomes globally chaotic. 

Coherent state representation — . In view of the large- 
ness of the spin, j ^ 1, we find it convenient to represent 
the theory in terms of coherent states 

\^)^j^^e^''-\j,j), (2) 

where z G C and \j,j) is a 'maximum- weight' cigcnstate 
of Jz, i.e. Jz\j,j) — The states \z) are unit nor- 

malized, {z\z) — 1. They entail the expectation values 
{z\Ja\z) — jla, where lx,y,z are the three components 
of a unit vector 1 = {sin 9 cos <j), sin 9 sin cf), cos 9)^ , whose 
angular orientation is defined through z = e**^ tan(6'/2). 
Each 1 2;) has minimum angular uncertainty, characterized 
by the solid angle 47r/(2j + 1) which defines a Planck cell 
on the unit sphere. The overcomplete set {\z)} provides 
a resolution of unity as 1 = / ^^^^.^ |z)(-z|- 
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Similarly, oscillator coherent states are defined as [H] 

|a) =e-""'/2e""'|0), (3) 

where a £ C and |0) is the vacuum, a|0) — 0. The states 
I a) assign a minimal uncertainty product to displacement 
and momentum such that these quantities are 'confined' 
to a single Planck cell located at a; oc (a + a*)^/H, p oc 
i{a — a*)\/~h within the classical phase space. The com- 
pleteness relation here reads 1 = ^ / da 

We next represent the system's time dependent den- 
sity operator p{t) in terms of a coherent state based 
quasiprobability density. Among the many possible 
choices the Glauber Q-function, 

Q{a,z)= {a,z\p\a,z) , (4) 

7r(i + zz* 

turns out to be the most convenient one by far, for 



our purposes. We here denote by \a,z) = \a) (g) \z) 
the overall coherent state for spin and oscillator. Us- 
ing the completeness relations given above one sees that 
Q is normalized as J dadz Q{a, z) — 1 and yields ex- 
pectation values of (anti-normal ordered) operators as 
(a"'at") = / dadza™a*"(9(a,z) . By its definition, Q 
exists and is non-negative, Q{a,z) > 0, for any density 
operator p. The latter property allows Q to converge to 
the classical phase space density as H ^ 0. As the most 
rewarding property we shall presently find that Q enables 
us to map the quantum evolution equation ihdtp = [H, p] 
onto a differential equation of Fokker-Planck type, opti- 
mally suited to describe the dissipation-free equilibration. 

Quantum dynamics — . Using identities such as 
l\z){z\ = ((1 - z^)d, + 2jj5+^)|z)(z| and a^\a){a\ = 
{da + a*) I a) (a I, it is straightforward to show that the 
evolution equation of the distribution function Q assumes 
the Fokker-Planck form 



C = ida ( ioa + g\/2j 



z + z* 



(5) 



+ id, (-cjo^ + ^(1 - z^){a + a*)] + c.c. , Cdm = 



/2j 



dad,{l - z^) + c.c. . 



h = ujo cos 9 + ujI + gV 81 cos ip sin 9 cos 



Most remarkably, derivatives terminate at second order. 
Less prudent choices of a quasiprobability might have 
brought the plague of higher derivatives [13] . 

To understand the meaning of the first order, or drift 
term, it is convenient to introduce canonically conjugate 
action angle variables {Ijip) as a = \/JI exp{iil)) for the 
oscillator and the canonical pair (cos 9, cj)) for the spin. It 
is then a straightforward (if tedious) matter to show that 
the drift operator becomes the classical Liouville opera- 
tor, £Q = {h, Q}, where the effective Hamilton function 



(6) 



obtains from the operator (jlj) by substituting J — >■ jl, 
~^ VJl G^vii'fp) and dividing out hj = L. The operator 
C thus describes the drift of the quasi-probability along 
the classical trajectories of the Hamiltonian flow. 

The distinctive feature of the classical dynamics is 
global chaos in the superradiant regime, g > Qc [I]. Re- 
ferring for an in-depth discussion to [H] , we here merely 
note that for g > gc and for excitations energies only 
slightly above the ground state bound Eq ~ — ?wo(— ) 
ergodic dynamics is observed. This is illustrated in the 
Poincare sections of Fig. [l] where the top/bottom row 
are for low {AE/\Eo\ ~ 0.20)/high {l^E/\Eo\ ~ 30) exci- 
tation. Both sections (e) reflect dominance of chaos and 
indicate equilibration. In particular, the low-excitation 
portrait (e) reveals chaos as dominant, even though the 



energy shell includes but a small part of the Bloch sphere. 

Quantum diffusion — . How does quantum mechanics 
interfere with the classical drift? Quantum mechanics en- 
ters the Fokker-Planck equation ([s]) through the diffusion 
operator £difi- A key feature of the diffusion operator is 
the absence of second order derivatives w.r.t. only oscil- 
lator (xi_2 = or spin (a;3^4 = cos6', ^), phase space 
variables. The 'diffusion matrix' D defining the opera- 
tor £diff = dxidx^Dij{x) thus possesses a block off- 
diagonal chiral structure, D = ( q ) ■ While there will 
be no need to spell out the explicit dependence of the 
chiral block d on the variables x = (/, f/;, cos 9, (/>), we em- 
phasize its smallness in The 'chiral' block structure 
of D entails a secular equation for the eigenvalues of the 
form A'' — trdd^ + detdc?''' = 0. The four eigenvalues 
of D thus come in two pairs i^/a, ±\/b where a, b are the 
eigenvalues of the non- negative 2x2 matrix dd^ . Each 
eigenvalue is associated with an eigenvector correspond- 
ing to diffusive spreading (-1-) or 'anti-diffusive' shrinking 
(— ). These four 'quantum' directions form a Cartesian 
frame wherein the four directions distinguished by the 
chaotic drift lie askew — one expansive (unstable), one 
contractive (stable) , and two neutral (along the flow and 
transverse to the energy shell). 

Since Q is guaranteed existence and positivity, quan- 
tum diffusion (along the eigenvectors of D with positive 
eigenvalues) must oppose the contraction along the clas- 
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FIG. 1. Poincare sections generated by monitoring the projection {lx,ly) of 1 in the southern hemisphere at fixed values of 
the phase ^. For each parameter value, g/gc, nine trajectories of different on-shell initial conditions are sampled. Upper row: 
energy Ae ~ 0.2|eoj above the ground state and values g/gc a) 0.2, b) 0.7, c) 0.9, d) 1.01, e) 1.5. Lower row: energy Ae ~ 20|eol. 



sically stable manifold, thus preventing the buildup of 
singular phase space structures. Deterministic contrac- 
tion and quantum diffusion balance at quantum scales 
1/a/J ^ V^- At shorter 'length' scales (which can be 
brought into play by initial states squeezed along the 
classically stable direction), the second order differen- 
tial diffusion operator becomes dominant (as can be seen 
by straightforvifard scaling estimates); it prevents further 
contraction and smoothes the distribution. By contrast, 
at the large scales generated in the classically unstable di- 
rection, the quantum (anti-) diffusive correction is of no 
significance. Rather, expansion continues uninhibited, 
accompanied by the chaotic 'folding' necessitated by the 
compactness of the energy shell. 

Toy model — . ft is instructive to illustrate the in- 
terplay of classical drift and quantum diffusion on an 
analytically solvable toy model, provided by the Fokker- 
Planck equation P = [dpp — dxX + D{dp — d1))P] the 
diffusion constant must be imagined small, < D <C 1 
and of quantum origin, D oc h. Anti-diffusive contrac- 
tion and deterministic stretching are carried by x ac- 
cording to va.rt{x) = (varo(a;) — £')e^* + D, where the 
variance is calculated using an averaging prescription 
{f{x,p)) ~ J dxdpf{x,p)P{x,p). We infer that the min- 
imum allowable 'quantum' scale of x is y/lD. If vai'tix) 
is of the order but larger than D, the characteristic scale 
for (the variance of) x will grow to ('macroscopic') order 
unity within the Ehrenfest time ^ ln(l/r'). On the other 
hand, diffusive expansion and deterministic contraction 
are carried by p as vart(p) = (varo(p) — D^e^^^ +D. The 
balancing scale for p is \/T). If the initial variance of p 
is much larger than D (say, of the 'macroscopic' order 
unity), the exponential decrease to the order of D takes 



a time of the order \n{l/D), again the 'Ehrenfest time'. 

Generalizing the toy model one may choose the axes 
of classical contraction and expansion skew to the prin- 
cipal axes of quantum (anti)diffusion. Positivity of all 
variances then restricts the relative orientations: clearly, 
the classically stable direction must not coincide with the 
axis of antidiffusion, or else the respective variance would 
go negative within a finite time. Referring for a detailed 
discussion of the ensuing correlations in the Dicke system 
to (14j , we here merely note that its guaranteed existence 
protects the Q-function from such type of alignment, save 
for exceptional regions in phase space. 

Thermalization — . Turning back to the full problem, 
we conclude that the interplay of deterministic expan- 
sion/contraction, quantum diffusion, and chaotic fold- 
ing (the latter of course absent in the toy model) will 
spread out any initial distribution over the compact en- 
ergy shell. Location, E, and width, AE of the shell 
are determined by the initial state, with ~ 
a minimal value for coherent state initial distributions. 
Given any fixed phase space resolution, Ax, and a char- 
acteristic width xq 3> Ax of the initial distribution, the 
quantum and a purely classical description of the flow, 
resp., both predict full coverage of the energy shell for 
time scales t > r ln(Aa;o/Aa;), where r is the Lyapunov 
time. However, important differences occur in the dy- 
namical buildup of the equilibrated state: while the clas- 
sical approach describes equilibration in terms of the for- 
mation of an infinitely filigree structure of alternating 
high and vanishing phase space density, — the result of 
continued stretching and folding of the initial distribu- 
tion — , the quantum distribution remains smooth (on 



scales 



J 



^1/2 



) all the way along. This hallmark of 
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quantum chaotic propagation should be experimentally 
visible by current spectroscopic techniques [15] • Finally, 
let us remark on the universality of the above picture 
with regard to variations in the initial conditions. Co- 
herent states, being isotropically 'supported' by a single 
Plack cell, come closest to the classical fiction of a sin- 
gle phase space point. They naturally qualify as initial 
states for the equilibration described above. Equally well 
suited are squeezed minimum uncertainty states, even if 
squeezed along the classically stable direction; the diffu- 
sion then smears the pertinent uncertainty to the char- 
acteristic quantum scale mentioned. Yet broader states 
succumb to equilibration even more willingly. 

Beyond the Dicke model — . While the relevance of 
'quantum diffusion' operators to the description of the 
long time dynamics in chaotic quantum systems has been 
noted before, previous work |l6i has added these contri- 
butions by hand. Our present analysis exemplifies how 
quantum diffusion emerges naturally. Evolution equa- 
tions with derivatives terminating at second order are 
not an exclusive privilege of the Dicke model. Rather, 
whenever a chaotic dynamical system has a Hamiltonian 
of the form of a second-order polynomial in the perti- 
nent observables and allows for a coherent-state-based Q- 
function, we expect a Fokker-Planck equation to govern 
the time evolution of Q and everything to go through in 
much the same way as above. Examples are autonomous 
SU{3) dynamics [T7] whose Hamiltonians contain terms 
of first and second order in the SU{3) generators. Cer- 
tain kicked systems qualify as well. Most notable among 
those is the kicked top [IBJ whose near classical quantum 
behavior has recently been observed experimentally ^15j . 
Genuine many-body systems also have Fokker-Planck 
equations for Q, provided the Hamiltonians are quartic 
in creation and annihilation operators but contain no an- 
tiresonant terms. Examples of much current interest are 
Bose- Hubbard systems [ini [20] • Finally, chaotic dynam- 
ics where the evolution equations of Q contains higher- 
order derivatives can allow for reasonable Fokker-Planck 
approximations, if good arguments are available for ne- 
glecting third and higher derivatives. We shall cover such 
extensions of our work in a separate publication |T3] . 

Summarizing, we have discussed the principles 
whereby the probability flux of a classically chaotic quan- 
tum systems covers its shell of conserved energy. The dy- 
namics crucially hinges on effectively diffusive quantum 
propagation interfering with the deterministic exponen- 
tial contraction inherent to classically chaotic flows. (In 
the expansive directions of the flow, quantum corrections 
are negligible.) The net effect of this competition is the 
smoothing of a flow which in a classical system would 
soon turn into a singular structure. We discussed these 
phenomena on the example of the Dicke model where the 
(Husimi representation of the) quantum flow equation 
assumes the particularly handy form of a Fokker-Planck 
equation. While the Fokker-Planck representation per- 



tains to other quantum systems, it is far from generic. 
In general, one will meet higher-order, differential equa- 
tions. However, power counting arguments |14j suggest 
that at the microscopic length scales ~ h}^^ where quan- 
tum diffusion becomes effective, higher terms in the ex- 
pansion will be sub-dominant. It thus stands to reason 
that the smoothing of the classical flow by mechanisms 
similar to those discussed above is a general effect. 
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